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Abstract 

In this paper, we have investigated the correlated electron 
emission of the nonsequence double ionization (NSDI) in an 
intense linearly polarized field. The theoretical model we em- 
ployed is the semiclassical rescattering model, the model atom 
we used is the helium. We find a significant correlation be- 
tween magnitude and direction of the momentum of two emis- 
sion electrons, and give a good explanation for this striking 
phenomenon by observing the classical collisional trajecto- 
ries. We argue that this correlation phenomenon is universal 
in NSDI process, as revealed by the recent experiment on the 
argon. 

PACS numbers: 32.80.Rm, 42. 50. Hz, 31.15.-p 



The excessive double ionization observed in Helium ex- 
periments |l]-|3| draws much attention to the multiple- 
electron dynamics in the laser-atom interaction. In 
these experiments the single ionization of He in a lin- 
early polarized field is accurately predicted by the sin- 
gle active electron (SAE) approximation B, well de- 
scribed by the Ammosov-Delone-Krainov (ADK) tun- 
nelling theory _|]. However, the case of double ionization 
is more complicated. In the regime of very high inten- 
sities (/ > 10 16 W/cm 2 ) where strong double ionization 
occurs, the double ionization keeps in good agreement 
with the sequential SAE models as that in the lower in- 
tensities regime(I < 10 14 W/cm 2 ). The double ionization 
deviates seriously from the sequential SAE model and 
shows a great enhancement in a "knee" regime [(0.8 — 
3.0) x 10 15 W/cm 2 ]. This surprising large yields of the 
double ionization obviously indicates that the sequential 
ionization is no longer the dominating process in this 
regime and the electron-electron correlation has to be 
taken into account. Intense efforts to model the two- 
electron process of the double ionization in a laser field 
have reproduced the main feature of the knee structure 
in the double ionization yield as a function of laser peak 
intensity and, moreover, yielded quantitative agreement 
with the experiments in some cases |^|.[| . 

The physical mechanism behind this nonsequential 
process is, however, still debatable. Both the "shake-off" 
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model and the "recollision" model are suggested to de- 
scribe the electron's correlation J_|,__,_]_| • However, none 
of the two nonsequence double ionization (NSDI) mech- 
anisms can completely explain the experimental obser- 
vations. For the "shake-off" model, it can not give the 
reason for the decrease in the double ionization yields 
as the polarization of the laser field departs from linear 
|p|-|rTf . In the "recollision" model, the returning electrons 
are known to have a maximum classical kinetic energy 
of ~ 3.2U P (Up = e 2 F 2 / 4m e ui 2 ) , so one can determine 
a minimum intensity required for the rescattering elec- 
tron to have enough energy to excite the inner electron. 
But the double ionization yields observed in experiments 
has no such an intensity threshold. In fact, the dou- 
ble ionization process is rather complicated and subtle, 
both of the two NSDI processes and the sequential ion- 
ization have contributions to it and may dominate in the 
different regimes. In another aspect, Becker and Faisal 
proposed a "correlated energy sharing" model to describe 
the NSDI processes and nuclei recoil experiment ^,0,1^ . 
The model is based on the so-called intense-field many- 
body S-matrix theory derived by a rearrangement of the 
usual S-matrix series and include time electron correla- 
tion and the rescattering mechanism. 

We have employed the semiclassical model to study the 
double ionization of helium in intense linearly polarized 
field (__,___]■ Our calculations reproduced the excessive 
double ionization and the photoelectron spectra observed 
experimentally both quantitatively and qualitively, and 
we argue that the classical collisional trajectories is the 
main source of the nonsequence double ionization of he- 
lium in the "knee" regime. 

Recently, the observation of the correlated electron 
emission in laser-induced double ionization of argon Jll| 
provided new insights into the NSDI process. These au- 
thors reported a strong correlation between the direction 
and the magnitude of the momenta of two electrons emit- 
ted from an argon atom: the momenta of the two emis- 
sion electrons tend to have the same magnitude and sign 
in the polarization direction. On the theoretical side, by 
solving the time-dependent Schrodinger equation for two 
electrons in three dimensions, Taylor et al [^6) gave that 
the most of double ionization probability flux tends to 
emerge to the same side of the ion. Similar conclusions 
have been drawn from one-dimensional model [jl7|| . 

In this paper, based on the 3D semiclassical rescatter- 
ing model developed recently [§_,___], we investigate the 
dynamical behavior of the correlated electrons in the dou- 
ble ionization process by analyzing their classical trajec- 
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tories. This investigation, as shown later, is very helpful 
to understand the physical mechanism behind the mo- 
mentum correlation. The model atom we use is the he- 
lium, however we argue that our discussions are available 
to the other multi-electron atoms, like the argon in recent 
experiment 15|. 

Firstly, we briefly present the semiclassical rescatter- 
ing model adopted in our calculations. The ionization of 
the first electron from bound state to the continues state 
is treated by the tunnelling ionization theory generalized 
by Delone et al. pq |. The subsequent evolution of the 
ionized electron and the bound electron in the combined 
Coulomb potential and the laser fields is described by a 
classical Newtonian equation. To emulate the evolution 
of the electron, a set of trajectories is launched with ini- 
tial conditions taken into from the wave function of the 
tunnelling electron. 

The evolution of the two electrons after the first elec- 
tron tunnelled are described by the classical equations 
(in atomic unit): 



d 2 r t 
dt 2 



-V(V£ + Vee) ~ F(t), i=l,2. 



(1) 



Here F(t) = F cos(ut)e z is the laser field. The indices 
i = l and 2 refer to the tunnel ionized and bound electron 
respectively. The potentials are 



K - v. 
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The initial condition of the tunnelled electron, under 
the S AE approximation of He + , is determined by a equa- 
tion including the effective potential given in Ref. Q 
and a generalized tunnelling formula developed by De- 
lone et al. [fill. In parabolic coordinates, the Schrodinger 
equation for a hydrogen-like atom in a uniform field e is 
written (in atomic unit), 



+ (-£ + i~ + TT + 1^)4 = 0, 
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in which I p i = —0.9 a.u. is the negative ionization po- 
tential of the outer electron. 

The above equation has the form of the one- 
dimensional Schrodinger equation with the potential 
U(rj) = -l/4i7 - I/877 2 - e-q/8 and the energy K = 

The evolution of the outer electron is traced by launch- 
ing a set of trajectories with different initial parameters 
to and vixo, where v± x q is the initial velocity perpendic- 
ular to the polarization of the electric field. The ini- 
tial position of the electron born at time to is given by 
x io = 2/io = , zio = —770/2 from the Eq. (||). The 
initial velocity is set to be v% y o = Vi z q = 0, Vx x q = v±o- 
Thus, the weight of each trajectory is evaluated by |18fi 



w(t ,v w ) = w(0)w(l), 



w(l) = 



y/2I pl V W 



ex P(-V /2/ pi w io/ e )> 
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and where w(0) is the tunnelling rate in the quasistatic 
approximation ]2(J. 

The initial state of the bounded electron is described 
by assuming that the electron is in the ground state of 
Hc+ with energy Ei = —2.0 a.u. and its initial distribu- 
tion is microcanonical distribution pijpf ■ 

In our calculation, the Eqs. (|]) are solved in a time 
interval between to and 15T by employing the standard 
Runge-Kuta algorithm. During the first ten optical cy- 
cles the electric field amplitude is constant, and then the 
field is switched off using a cos 2 envelope during three 
cycles, and during the last two optical cycles the elec- 
trons is free from the electric field. The wavelength is 
A = 780 nm, which is so chosen to match the experiment 
^Hl, an d the intensity of laser is I = 1 x l0 15 W/cm 2 
which is right at the 'knee' region of the double ioniza- 
tion of helium. 10 6 or more initial points are randomly 
distributed in the parameter plane —ir/2 < uto < 3ir/2, 
vixo > for the outer electron and in the microcanoni- 
cal distribution for the inner electron. The distribution 
for the ionization electron can be obtained by making 
statistics on an ensemble of classical trajectories weighed 
by the (^). The results have been tested for numerical 
convergence by increasing the number of trajectories. 

Figure 1 shows the momentum correlation between the 
two emission electrons in the double ionization of the 
present calculations. The horizontal axis shows the mo- 
mentum component of the first electron in the direction 
of polarization (Pi z ) and the vertical axis the same mo- 
mentum component of the second electron {Piz)- This 
figure shows a strong correlation between the momenta 
of the two electrons. There is a clear maximum for both 
electrons being emitted with the same momentum com- 
ponent in the direction of polarization axis of about 2.7 
a.u., and emission to opposite half planes is strongly sup- 
pressed, i.e., both two electrons tend to fly to same side of 
ion in the direction of polarization. This phenomena has 
been observed in the 'knee' region for argon JT|. On the 
other hand, from Fig. 1, we see that the maximum mo- 
mentum of both electron is about 4.5 a. u., which is consist 
with the electron-ion coincidence experiment observation 
of helium p3[ | , in which the maximum energy of emission 
electron in NSDI process is 4U p , since the perpendicular 
component of momentum is small, the maximum mo- 
mentum component in the polarization direction can be 
approximate obtained as P zmax = \/8Up ~ 4 a.u.. 

An useful alternative perspective on Fig. 2 is obtained 
by rotating the distribution by 45° . Then we can get two 
new distributions. In Fig. 2(a), we show the distribution 
of the sum momentum P + = (Pi z +P2 Z )', In Fig. 2(b), we 
show the distribution of the difference momentum P~ = 
(Piz — Piz)- Owing to momentum conservation, the sum 
momentum of emission electrons is equal and opposite 
to the He 2+ recoil ion momentum P2|jl4| l, so its distri- 
bution, as we know, shows a characteristic double-hump 
structure with a central minima. A vast amount of litera- 
ture has devoted to discuss the double-hump structure of 
momenta of recoil ions f2^,[l4|] , here we only pay attention 
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to the correlation between the two emission electrons. 

Comparing the Fig. 2 (a) with Fig. 2(b), one finds that 
the sum momentum at the peak is about 4.3 a. it., al- 
most above the maximum momentum of each electron, 
and the maximum sum momentum is almost 8.7 a.u., 
about twice of the maximum momentum of one electron; 
furthermore the distribution width of difference moment 
is much smaller than the one of the sum momentum. 
These features indicate the momenta parallel to polar- 
ization axis of the two emission electrons likely have the 
same direction. On the other hand, the peak of the dis- 
tribution of difference momentum is at zero, so the mo- 
menta of the two emission electron tend to have the equal 
magnitude. 

To study the origin of the correlation of electrons emis- 
sion, we show two trajectories of the electrons in Fig. 3. 
Fig. 3(a) shows a typical trajectory of which the dif- 
ference momentum in polarization direction is small and 
fig. 3(b) shows a typical trajectory of which the difference 
momentum in polarization direction is nearly the maxi- 
mum. From these trajectories, it can be concluded that 
after the second electron ionized the momentum trans- 
fer received from the laser field is almost identical for 
both electron. Because the electron can obtain very high 
energy from the laser, the momenta of the two thus be 
accelerated to the same direction. On the other hand, 
this part of the acceleration only adds to the sum of the 
momenta of the two electrons, but does not change the 
difference of the momenta. Therefore the difference of the 
momenta is only determined by the ionization process of 
the second electron. 

As we have discussed before, the second electron is 
mostly ionized by the "recollision" mechanism: the sec- 
ond electron is ionized by a collision with the tunnelled 
electron Q. Because the collision between the two elec- 
trons is almost instantaneous so that the energy is con- 
served approximately when the collision happens. The 
total energy of the two electron can be expressed by 

E = E T+ {^^)\{^^ 2 (6) 

where Et includes the potentials and the kinetic energy 
of the perpendicular parts. To the best of our knowl- 
edge, so far the energy distribution of the returning elec- 
tron has not been reported, but we believe that it ex- 
hibits a peak at zero energy and decreases as the energy 
increasing, on the other hand, only the returning elec- 
tron which energy is above the ionization potential of the 
bound electron I P 2 can cause the second electron ionized, 
so the total energy of the two electrons soon after the col- 
lision occurred is small and its distribution peaks at zero. 
Therefore the P~ is more likely zero, i.e., the momen- 
tum components in the polarization direction of the two 
emission electrons are likely equal to each other. The re- 
turning electrons are known to have a maximum classical 
kinetic energy of ~ 3.2U P , therefore the maximum differ- 
ence momentum must satisfy ( Plz ~ 2 P ' 2 * ) ~ 3.2U P — I P 2, 



so \P~ | max = 2^/2,.2U p - I p2 ~ 4.2 a.u.. 

On the other hand, since the total energy of the two 
electrons soon after the collision occurring is small, the 
total final energy of the two emission electrons mostly 
reflect the energy transfer received from the field. The 
field acceleration make the sum momentum increase. The 
maximum energy for an electron in the double ionization 
process is about 4U P from Eq. (||), assuming the two 
emission electron have the same energy, we can obtain 
the maximum sum momentum |P + | max — 8 a.u.. This 
indicates that the joint acceleration of the electrons in 
the laser field clearly dominates over the influence of the 
electron repulsion, and both electrons ionized are driven 
by the laser to the same side. 

In conclusion, we have investigated the momentum 
correlation between magnitude and direction of the two 
emission electrons in nonsequence double ionization. The 
numerical results on the helium show a significant corre- 
lation on the momentum of the two electrons: the emis- 
sion electrons tend to have the same momentum com- 
ponent in the polarization direction. The phenomena 
can be directly comprehended from the classical colli- 
sional trajectories. These discussions suggest that the 
correlated electrons emission in double ionization pro- 
cess is a semiclassical process. We also evaluated the 
width of distribution of the sum and difference momen- 
tum. Because the difference momentum is only deter- 
mined by the ionization process, so it is important to ver- 
ify the dominating process in the "knee" regime. Based 
on the rescattering model, we argue that the maximum 
difference momentum of the two emission electrons is 
\P~ | max = 2^/3. 2U V — I V 2- The predictions coincide with 
the argon experiment Ea . We hope our discussions will 
stimulate the experimental works on the helium in the 
direction. 
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FIGURES CAPTION: 

Fig. 1. Momentum correlation between the two emit- 
ted electrons given by present calculations. 

Fig. 2. The distribution of momentum parallel to the 
polarization axis: (a) the sum momentum of the two 
emission electrons, (b) the difference momentum of the 
two emission electrons. 

Fig. 3. Two typical trajectories of the " rccollision" 
mechanism: (a) the case that the difference momentum 
is small, (b) the case that the difference momentum is 
nearly the maximum. 
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